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Abstract. One of the most interesting and current phenomenological extensions of General 
Relativity is the so-called f(R) class of theories; a natural generalization of this includes an 
explicit non-minimal coupling between matter and curvature. The purpose of this work is 
to present a unified view of the applicability of the latter to various contexts, ranging from 
astrophysical matter distributions to a cosmological setting. Various results are discussed, 
including the impact of this non-minimal coupling in the choice of Lagrangian density, a 
mechanism to mimic galactic dark matter and a Cosmological Constant at a astrophysi- 
cal scale, the possibility of accounting for the accelerated expansion of the Universe and 
modifications to post-inflationary reheating. The equivalence between a model exhibiting a 
non-minimal coupling and multi-scalar-theories is also discussed. 
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1. Introduction 

Cosmology currently faces two outstanding 
challenges: the origin of dark energy and dark 
matter Research on these dark components 
of the Universe may be roughly divided into 
two fields: one assumes that Einstein's General 
Relativity (GR) is valid at the relevant en- 
ergy scales, thus focusing on the search for 
the missing matter or energy contribution; the 
other posits instead that dark matter and energy 
merely reflect deviations of gravity from GR. 

When addressing dark matter, the first ap- 
proach relies fundamentally on the character- 
ization of additional matter species, arising 
mostly from extensions to the Standard Model 
of particles and fields, collectively dubbed as 
weak-interacting massive particles (WIMPS): 
these include e.g. neutralinos, which arise nat- 
urally from supersymmetry, or the axion, de- 



rived from an elegant solution to the strong CP 
problem (Feng 2010). 

Similarly, the nature of dark energy usu- 
ally involves the addition of a new scalar field, 
which slow-rolls down an adequate potential 
— quintessence (Copeland et al. 2006). A pu- 
tative unification of dark matter and energy has 
also been suggested, by resorting to a scalar 
field (Bertolami & Rosenfeld 2008) or an ex- 
otic equation of state (EOS), such as the gen- 
eralized Chaplygin gas (Kamenshchik et al. 
2001; Bilic et al. 2002; Bento et al. 2002). 

The second approach follows a quite dif- 
ferent route, by assuming that GR itself is at 
fault: instead of dark matter or energy, one may 
consider that dark gravity is at play. Thus, sev- 
eral attempts to generalize GR have surfaced: 
these are usually phenomenological in nature, 
positing putative low-energy modifications of 
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Einstein's theory that should be derived from 
a yet unknown high-energy, fundamental the- 
ory of gravitation. In a cosmological context, 
modifications of the Friedmann equation to in- 
clude higher order terms in the energy density 
p have been proposed (Maartens 2004; Freese 
& Lewis 2002; Carroll et al. 2004), as well as 
considerations of the impact of a van der Waals 
EOS for matter (Capozziello et al. 2002). 

Another widely discussed candidate is the 
Modified Newtonian Dynamics (MOND), a 
non-relativistic modification of the Poisson 
equation, which becomes V[yu(V<l)/ao)V<l'] - 
4nGp, where <1> is the gravitational potential, 
p the density and ^(x) a function satisfying 
fi(x) ~ X for X « 1 and fi{x) ~ 1 for x » 1 . 
By introducing the typical acceleration ao ~ 
10"'° m/s'^, MOND allows a simple derivation 
of the flattening of the rotation curves of galax- 
ies, thus accounting for the missing "dark" 
matter challenge (Milgrom 1983). 

However, MOND has been shown to suf- 
fer from one outstanding issue, namely the in- 
ability to explain the gravitational profile of the 
so-called Bullet cluster (Clowe et al. 2006) — 
although a proposal resorting to heavy neutri- 
nos with a mass close to ~ 2 eV has been 
put forward (Angus et al. 2006, 2007). On 
the theoretical side, the current relativistic for- 
mulation of MOND, dubbed Tensorial- Vector- 
Scalar (TeVeS) theory, resorts to two scalar 
(one dynamic and one auxiliary) and one vec- 
tor field, together with a physical metric gf^y 
used to couple with matter (Bekenstein 2004): 
this formulation may scarcely be considered an 
elegant one, and one is tempted to call these ex- 
tra fields as dark matter in another guise. 

Changes to the action functional provide 
a simple pathway for implementing exten- 
sions to GR: a straightforward approach re- 
places the linear scalar curvature R term in 
the Einstein-Hilbert action with a function 
/(/?) (Capozziello et al. 2005; Allemandi et al. 
2004) (see (Felice & Tsujikawa 2010) for a re- 
view); more evolved dependences may also be 
explored, as considered in Gauss-bonnet mod- 
els (Lovelock 1971). 

This class of f(R) theories enjoys con- 
siderable success in several fronts, including 
the puzzle of the missing "dark" matter in 



galaxies and clusters (Capozziello et al. 2007, 
2009a), as well as the nature of "dark" en- 
ergy (Capozziello et al. 2003); the early pe- 
riod of rapid expansion of the universe is de- 
scribed by the Starobinsky inflationary model 
/(/?) =R+aR^ (Starobinsky 1980); local Solar 
system impact and the related parameterized 
post-Newtonian (PPN) metric formalism have 
also been studied (Capozziello et al. 2009b). 

The mathematical equivalence between 
f(R) and scalar-tensor theories has proven very 
profitable, given the possibility of comparing 
results between the two types of model (Felice 
& Tsujikawa 2010). Further insight is gained 
from considering the so-called Palatini formu- 
lation, where both the metric and the affine 
connection are taken as independent variables 
(Sotiriou et al. 2006) — an alternate to the 
usual metric affine connection approach, where 
the affine connection is taken a priori as de- 
pending on the metric. 

Given the proficuous results arising from 
f(R) theories, one is naturally tempted to fur- 
ther generalize this model — in the process 
extending its explanatory capability. Thus, an- 
other interesting possibility has been subject 
to scrutiny: not only that the curvature is non- 
trivial in the Einstein-Hilbert Lagrangian, but 
also that the coupling between matter and ge- 
ometry is non-minimal (Bertolami et al. 2007), 
as written below: 

5 = J lKR+f2(R)£] ^d'x . (1) 

where k = c'^/167tG. 

The purpose of this work is to present an 
unified view of several results that arise from 
this framework in a wide range of scenar- 
ios (Bertolami & Paramos 2008, 2010, 2011; 
Bertolami et al. 2008, 2010, 201 1) ; these were 
developed in collaboration with O. Bertolami, 
F. S. N. Lobo and P. Frazao. 

2. The model 

Variation with respect to the action Eq. (1) 
yields the modified Einstein field equations, 

2 (kFi + F2£) G^y = 2A^v {kF^ + F^L) (2) 
-Vk{FxR - /i) + F^mUy + h^y^y - 
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where one defines A^v = ^fi^v - guvO and 
FiiR) = fliR). As expected, GR is recovered 
by setting fi(R) - R and f2(R) = 1. By taking 
the trace of the above, one gets 



KiFiR - 2/0 + F2£R = 
^f2T-3DiKFi+F2£) , 



(3) 



where T is the trace of the energy-momentum 
tensor. The usual f(R) theories are recovered 
by setting f2iR) = 1, so that Eq. (3) yields an 
algebraic relation for R = R(T). However, if 
one considers a non-minimal couphng /2(/?) i= 
1, the above becomes a difl'erential equation: 
in particular, this enables the possibility that, 
given a rapidly varying F2JI, the above may 
lead to high curvatures even if T is low. 

2.1. Non-conservation of energy 

The Bianchi identities may be used to derive 
the non-(co variant) conservation of the energy- 
momentum tensor, 



h 



(4) 



In the absence of a non-minimal coupling, 
fliR) = 1. one recovers the covariant conser- 
vation of the energy-momentum tensor. 

Eq. (4) may be rewritten as an extra force 
imparted on test particles, so that its trajectory 
will deviate from geodesical motion. Thus, the 
Equivalence Principle may be broken if the 
r.h.s. of the last equation varies significantly for 
difi'erent matter distributions. 

In Jordan-Brans-Dicke theories, the curva- 
ture appears coupled to a scalar field (Damour 
& Esposito-Farese 1992), leading to expres- 
sions similar to Eq. (4), due to the en- 
ergy exchange between matter and the latter. 
However, this non-conservation may be trans- 
formed away by a suitable conformal transfor- 
mation to the Einstein frame, where the curva- 
ture appears uncoupled (Faraoni et al. 1999). 
By the same token, one could assume that the 
non-conservation law Eq. (4) is only an artifact 
of the particular frame adopted, not a physi- 
cally significant feature of the model. 



However, performing the conformal trans- 
formation gf^y g^y = f2gijv, one finds that en- 
ergy conservation is attained, Vftff^ - 0, only 
if f"" = f-^T"" and 2£ = T. Unless this re- 
lation holds for all matter species, energy con- 
servation is indeed broken. 

Thus, it suffices to use one example that 
breaks this relation: taking the case of a per- 
fect fluid, with a Lagrangian density Jl = -p 
(Bertolami et al. 2008), then one obtains p = 
-p/3 — signaling a perfect fluid with nega- 
tive pressure (but different from a cosmologi- 
cal constant, where one has - -Pa); if one 
instead opts for the Lagrangian density £, = p 
(Sotiriou & Faraoni 2008), then p = p is ob- 
tained — the EOS for ultra-stiflF matter. 

The discussion of which form for £, is cor- 
rect in the presence of a non-minimal coupling 
is addressed in the following section; never- 
theless, it is clear that one cannot recover en- 
ergy conservation for a perfect ffuid — and, 
by extension, for all matter species: this non- 
conservation law is indeed a fundamental prop- 
erty of the model under scrutiny. 



3. Equivalence with 

multi-scalar-tensor models 

Through a suitable conformal transformation, 
the usual f(R) theories can be rewritten as GR 
with an added scalar field contribution, which 
is dynamically identified with the curvature, 
(f> = R. Similarly, the discussed non-minimally 
coupled model Eq. (1) can be recast as a multi- 
scalar field theory, with two scalar fields, al- 
beit a "physical" metric remains in the mat- 
ter Lagrangian (Bertolami & Paramos 2008). 
Indeed, by performing a conformal transfor- 
mation gf^y -» gf,y = exp[(2/ y/3)(p^]gijy, Aic 
equivalent action is obtained. 



f2{ip^)e-Ty£(g^^,;^) (5) 



+2K{R-2g^(Tijip\^>fily-AU) 
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where X denotes all matter fields, (p^ and (p^ are 
scalar fields related to the scalar curvature and 
matter Lagrangian density through 



1 V3, 
<P = — log 



Fi(R) + F2{R)£ 



2k 



and (fi^ = R; CTij is the field-metric 



(Tii = 



1 1 

-1 



and the potential is given by 



I 7 1 

U{ip\(p') = -exp 



2k 



■ exp 



3 

2V3 J 



(6) 



(7) 



(8) 



Notice that any use of the metric in this 
Lagrangian density (contractions in kinetic 

terms, etc.) is done not with f^y, but with the 

physical metric g^y = exp[-(2/ V3)v?']^jl<v 

4. Lagrangian density of a perfect 
fluid 

The argument of the preceding section serves 
to show another striking feature of non- 
minimally coupled models: the Lagrangian 
density of matter appears explicitly in the field 
equations. In particular, one focuses on the 
case of a perfect fluid, given its ubiquity as 
a useful description for standard matter. Its 
energy-momentum tensor is the familiar form 



(P + p)UftUy + Pgfjy , 



(9) 



where p is the energy density, p the pressure 
and the four- velocity (with m^m'' = -1); its 
trace is J = 3p - p. The associated Lagrangian 
density is less discussed, precisely because it 
is absent from the field equations of GR: in 
applications of GR, the choice of Lagrangian 
density is mostly irrelevant, while it is of the 
utmost importance when considering a non- 
minimal coupUng (Bertolami et al. 2008). 

The identification Jl - p was first ad- 
vanced in (Sehger & Whitham 1968), with 
a relativistic generalization in (Schutz 1970). 



Much later, (Brown 1993) showed that this 
choice is equivalent to X = -p, complemented 
by a suitable set of thermodynamical potentials 
and Lagrangian multipliers ^ - ip^^ + s9^^ + 
PAO^fx', these enable the relativistic thermody- 
namical relations via a current term J^(p^ in the 
action, where is the vector density, i.e. the 
flux vector of the particle number density. In 
(Hawking & Elhs 1973), an isentropic perfect 
fluid is described via X = -p. 

The equivalence between occurs on-shell, 
by substituting the field equations derived from 
the matter action 



(10) 



Sm^ j {-^Fgp + J^(lf)(fx , 

back into the action functional and reads the 
resulting on-shell Lagrangian density Jl\ = p. 
Similarly, the action may be rewritten so that 
the on-shell Lagrangian density reads X2 - na, 
where n - \J\/ is the particle number den- 
sity and a{n, T) - p{n)ln - sT, where s is 
the entropy per particle and T is the temper- 
ature. Lastly, one may remove the current term 
■^/'<^/ from Eq. (10), thus obtaining the 
on-shell Lagrangian density X3 = -p through 
the addition of adequate surface terms. 

Indeed, (Seliger & Whitham 1968) did not 
argue that Hi - p was the bare (as opposed 
to on-shell) Lagrangian density for a perfect 
fluid, but resorted to a much more evolved ac- 
tion functional: this simple identification with 
the pressure also arose only after the resulting 
field equations where substituted into it. 

With the above in mind, one now ascertains 
how should this procedure be generaUzed to a 
non-minimally coupled scenario. This should 
afi'ect the terms in Eq. (10) that are minimally 
coupled, so that the matter action becomes 

S'm^ J{- ^F8f2(R)p + J't^,) d'x , (11) 

while the current term remains uncoupled 
(aside from the use of the metric to contract 
indexes). Varying the action with respect to the 
potentials included in 0^, one obtains 



n OS 



(12) 
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where jx is the chemical potential and 6 is a 
scalar field (included in 0^) whose equation 
of motion imposes the entropy exchange con- 
straint (sJ^)^fi - 0. Thus, the non-minimal cou- 
pling of curvature to matter is reflected in both 
the velocity and the temperature identification. 

By substituting the modified equations of 
motion into action (11), an on-shell Lagrangian 
density JLi= p may be read, as in GR. The ad- 
dition of surface integrals also yields the dis- 
cussed JL2 - -p and X3 = -na. 

Albeit the action may adopt distinct on- 
shell forms, this does not translate into a direct 
equivalence between them: only the original 
bare Lagrangian density Xo should be inserted 
into the field equations Eq. (2). However, 
this bare Xo should not appear into the non- 
conservation law: indeed, when deriving Eq. 
(4), the current term (which is not coupled with 
the metric) is dropped and one is indeed left 
with X2 = -p- 

Notwithstanding the above, one can adopt 
a simpler stance regarding the choice of the 
Lagrangian density £. - -p instead of £. - p: 
if a dust distribution is to be considered — 
that is, a perfect fluid with negligible pressure 
and corresponding energy-momentum tensor 
r^v = pg^y — , it appears unnatural to take a 
vanishing quantity as Lagrangian density. 

Ending this section, one remarks that if a 
matter form is actually described by two inde- 
pendent bare Lagrangian densities leading to 
diflerent dynamical behaviour of Eq. (2) and 
(4), then only confrontation with observation 
can ascertain which is the correct description. 

5. Dark matter mimicking 

One now reviews a mechanism that mimics the 
missing dark matter revealed through the flat- 
tening of galaxy rotation curves (Bertolami & 
Paramos 2010). In order to isolate the effect of 
the non-minimal coupling, one sets fi{R) = 1 
and assumes a power-law 



(13) 



Since the dark matter contribution is dominant 
at large distances i.e. low curvatures, a negative 
exponent n is expected. 



Inserting this into the modified field Eq. 
(2), together with the Lagrangian density £, = 
-p for a dust distribution with p = 0, yields 



P_ 
kR 



1 



(14) 



Rn 





n -I 


[ru 


- 

1 kR 



At large distances, normal matter is subdomi- 
nant and the trace Eq. (3) reads 



p_ 

kR 



■(15) 



Inspection shows that an exact solution is ob- 
tained if the last term vanishes. 



R - Rn 



ll/(l-n) 



(1 -2n)- 



Pn 



(16) 



defining the characteristic density p„ - 2kR„. 

A more evolved study of Eq. (15) shows 
that the most general solution oscillates around 
the one above: this makes the gradient term ac- 
tually dominate Eq. (15), and also allows for 
a perturbative f2(R) 1. As numerical re- 
sults show, one may disregard these oscilla- 
tions and simply consider Eq. (16) (Bertolami 
& Paramos 2010). A perturbative non-minimal 
coupling is paramount: it makes the mimicking 
mechanism satisfy both that the weak, strong, 
null and dominant energy conditions; grants 
immunity against Dolgov-Kawasaki instabili- 
ties (Dolgov & Kawasaki 2003); and makes the 
extra force arising from Eq. (4) very small. 

Instead of solving Eq. (14), obtaining a 
modified gravitational potential and then read- 
ing the mimicked dark matter contribution, 
one interprets the additional curvature obtained 
from the non-minimal couphng as due to the 
latter's density profile. 



Pdm = 2kR = Pn 



(l-2n)^ 

Pn 



l/(l-n) 



(17) 



Thus, one obtains a direct translation between a 
visible matter profile p and the mimicked dark 
matter distribution pj^, with n related to the 
large radius behaviour of both contributions. 

The mimicked dark matter may be further 
characterized by inserting Eq. (16) into Eq. 
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(14), and reading the obtained terms as due 
to a corresponding energy-momentum tensor: 
this allows one to obtain the EOS parameter 
for dark matter, w = PdmlPdm - n/il - n): for 
a negative exponent n, one has w < 0, hinting 
at at possible cosmological role in enabling an 
accelerated expansion of the Universe, as will 
be discussed in the following section. 

In order to fit the galaxy rotation curves 
of several galaxies (NGC 2434, 5846, 6703, 
7145, 7192, 7507 and 7626, which are almost 
spherical type EO galaxies with well deter- 
mined rotation curves (Gerhard et al. 2000)), 
one resorts to the Hernquist profile for visible 
matter (Hernquist 1990), which behaves as 



p(.)~^(i.^r . 



(18) 



and the Navarro-Frenk- White (Navarro et al. 
1995) and isothermal sphere profiles for dark 
matter: the latter enables a perfectly flat rota- 
tion curve, while the former is favored by nu- 
merical simulations: 

Pisir) ~ r-^ , PNFwir) ~ - (l + -) ' .(19) 

with a signalling the transition between in- 
ner and outer slope profiles. Given the outer 
slope behaviours p(r) ~ r p/s(r) ~ and 
PNFw(.r) ~ r^^, the relation pdm(r) °^ p(r)'^^'""' 
yields the exponents rijs = -1 and n^Fw = 
-1/3 that translate the Hernquist profile into 
the two above. Thus, one is led to consider the 
composite non-minimal coupling 



f2{R) = 1 + ^ 



R 



1/3 



R 



(20) 



The values of the model parameters R-\ and 
R-\/3 are obtained from fitting the result of 
the numerical integration of Eq. (14) to the 
available galaxy rotation curves: by doing so, 
one concludes that R^i ~ 1/(16.8 Gpcf and 
/?_i/3 ~ 1/(1.45 X 10^ Gpcf-. Fig. 1 shows the 
fit to NGC 5846, one of the galaxies addressed 
in (Bertolami & Paramos 2010). 

Allowing for individual fits of these pa- 
rameters to each of the considered galaxies 
shows that there is some dispersion around 
these orders of magnitude, which could be due 




Fig. 1. Observed rotation curve (dashed full), 
decomposed into visible (dotted) and dark mat- 
ter (dashed grey) contributions (Gerhard et al. 
2000), with mimicked dark matter profile (full 
grey) and resulting full rotation curve (full). 



to deviations from sphericity, poor choice of 
visible or dark matter density profiles, unac- 
counted effect of a non-trivial f\ (R) + R term 
or the non-minimal coupling with the electro- 
magnetic sector (thus enabling a dependence 
on luminosity) etc.. Nevertheless, the quality of 
the obtained fits and the elegance of the mim- 
icking mechanism support its ability to account 
for the missing dark matter puzzle. 



6. Accelerated expansion of the 
Universe 

This section discusses the possible use of a 
non-minimal coupling to describe the current 
phase of accelerated expansion of the Universe 
(Bertolami et al. 2010). One focuses on a con- 
stant deceleration parameter q = -aa/(a)^ = 
1/jS - 1, which translates into a power- law ex- 
pansion with the scale factor evolving as a{t) = 
aoit/tof, with to = 13.73 Gy. 

This hints at the use of a power-law for the 
non-minimal coupling also, as considered in 
the previous section. Thus, one takes /2(/?) = 
1 -I- (R/RnT, with the exponent n assumed to be 
negative — so that the accelerated expansion 
phase appears at late times, when 

This accelerated expansion is obtained in a 
rather straightforward fashion: assuming a flat 
Friedmann-Robertson-Walker metric with line 
element ds^ = -dt^ + a^dV^ and that matter is 
described by a perfect fluid with density p and 
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pressure p, one has Too = p and T^r = a^p. 
One then uses X = -p and Eq. (4) to ascer- 
tain that energy is covariantly conserved in a 
cosmological context, 



fi ^ ' 
p + 3Hp^Q^ pit) = po (y) 



(21) 



where po = SlmPcrit, with ~ 0.3 the rela- 
tive matter density and pcru ~ 10"^^ kg/rn^ the 
critical density (Spergel et al. 2007). 

One uses this result together with Eq. (2) to 
compute the modified Friedmann equation 



9 1 

= t(P + P+Pc + Pc) 

OK 



(22) 



where the additional density Pc and pressure Pc 
terms are introduced, 



-6pojS 



Pc = -2pon 



1-2/3 + n(5/3 + 2n - 3) 



,(23) 



(tr(^r [6/^(2/^-1)]- 
2 + 4rP-- P(2 + 3/3) + «(8/3 - 6) 
(tf(xf [6^2^-1)]- 



defining t„ = R^^^^ and taking the weak regime 
F2P « K (Bertolami et al. 2010). 

Given the Hubble parameter H(t) = a/a - 
/3/ 1, the l.h.s. of the Friedmann Eq. (22) falls as 
r^, so that comparing with the above gives a 
relation between the exponents jS and n: 




Fig. 2. Evolution of the deceleration parameter 
q(z) for n - -4, t2 - to /A (full) and n = -10, 
f2 = to/2 (dashed); Icr, 2cr and 3cr allowed 
regions are shaded, white hne gives best fit 
(Gong & Wang 2007). 



This value for the exponent n is very dis- 
tinct from the n/s = -1 and vimfw - -1/3 
scenarios considered in the previous section. 
These have no cosmological impact on the cur- 
rent scenario, given the broad difference be- 
tween the relevant timescale ?o and those ob- 
tained from f_i = RZ\^ and f_i/3 = RZ^^^l- Its 
value n = -10 may also be regarded as some- 
what unnatural, but is highly dependent on the 
q(z) profile: Fig. 2 shows that n - -4, t2 - to/A 
also yields a curve within the allowed 3cr re- 
gion. Despite this caveat, one concludes that a 
non-minimal coupling may be used to describe 
the accelerated expansion of the Universe. 



3fi + 2n = 2^/3= -(l-n) . 



(24) 



Thus, one concludes that any negative expo- 
nent n will yield an accelerating Universe with 
P, q> Q. One may also compute the decelera- 
tion parameter and equation of state parameter 
for the non-minimally coupled contribution. 



q = -\ + 



2(1 -«) 



1 



(25) 



The latter is the same form found in the previ- 
ous section. By fitting the numerical solution of 
Eq. (22) to the evolution profile of q{z) (Gong 
& Wang 2007), one obtains the best fit values 
n = -10 and ?_io = to/2, depicted in Fig. 2. 



7. Matter solution with constant 
curvature 

The main result of the previous section is em- 
bodied by Eq. (25), which directly translates 
the exponent n of a power-law coupling be- 
tween curvature and matter to the deceleration 
parameter q. However, it is trivial to check that 
one cannot produce a de Sitter Universe, char- 
acterized by <7 = -1, as this would require 
n — > 00 — so that a Cosmological Constant 
(CC) A with EOS pa = -pA cannot be recov- 
ered. One could argue that this hints at the need 
to consider a different form for f2{R) — per- 
haps an exponential, given the behaviour of the 
scale factor a{t) = uq exp{Hot). 
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However, as was shown in (Bertolami et al. 
2010), it turns out that a de Sitter phase is im- 
possible to obtain in a cosmological context, 
fundamentally because the curvature is given 
by the CC, R - 4A, and is thus constant: 
since the density still evolves as p{t) oc ~ 
exp(-3Hot), the terms arising from the non- 
minimal coupling in the modified Friedmann 
cannot be interpreted as an additional density 
and pressure due to a CC, pcc = -pcc- One 
may simplistically state that, since the f2(R)X. 
term in Eq. (1) cannot be made constant, a CC 
is unattainable. 

With the above in mind, this section briefly 
discusses the possibility that a CC is not di- 
rectly generated in a cosmological context, 
but arises due to an averaging of several 
constant curvature "bubbles" (Bertolami & 
Paramos 2011). This does not aim at solv- 
ing the Cosmological Constant (CC) problem: 
the ~ 10^^" discrepancy between its expected 
value and that needed to account for the cur- 
rent accelerated expansion, but draws a mech- 
anism to produce the latter (for a review, see 
Bertolami 2009). 

Notice that a constant curvature R = Rq 
implies that all terms involving /2 and its 
derivatives become constants; thus, the actual 
form of /2 is shown to be only constrained by 
^2 (4 A) A «c /2(4A). In particular, the linear or 
power-law non-minimal couplings considered 
in this work are compatible with this relation. 

One uses a static Birkhoff metric with line 
element Js^ = -e^''^''W+e^''^''^dr^+dn^; mat- 
ter is assumed to be described by a perfect fluid 
with constant EOS parameter o) = pjp. The 
Einstein field Eq. (2) with f\{R) = R thus leads 
to the diff'erential equation. 



''=f'"--f^'^WFFT2h^'^'- (26) 

r KKRorY F2 



_J_I l+OJ f(f2'^ 



2,r^\2{KRa)^j \F2 



1 



6(/s:/?o)V 



(p'f\--RQ + 



l+OJ f2\ 

2 F2] ' 



where K is an integration constant. Imposing 
continuity with a Schwarzschild-de Sitter met- 




Fig.3. Rescaled density profile arising from 
Rq = 4A. 

ric at the boundary of the matter distribution 
r = Vs, one obtains 



K 




(27) 



The boundary radius Tj is related to the total 
mass of the matter distribution via 



/— 11.4Mpc , 

Pc 



(28) 



where pn = 1 H atom/cm^ = 1 .660 X 
10"^^ kg/rrp. The (rescaled) numerical solu- 
tion to Eq. (26) when Rq = 4 A (i.e. the local 
curvature yields a CC) is shown in Fig. 3. The 
typical distance between galaxies is of the or- 
der of r, ~ IMpc (e.g. Andromeda is 778 Kpc 
away), so that pc ~ 130 H atom/cm-' — close 
to the lower bound of HII star-forming re- 
gions. The intersteUar medium attains densi- 
ties up to 10^ H atom/cm^, allowing for a ra- 
dius of the constant curvature bubble down 
to ~ 100 Kpc. Thus, one concludes that it is 
possible to obtain a smooth constant curvature 
matter solution R^ = 4A compatible with the 
expected values for realistic density profiles. 

7.1. Averaging meciianism 

Having obtained a feasible description of one 
constant curvature "bubble", one now briefly 
discusses how several of these may average 
and produce a CC at cosmological scales. One 
first argues that a universal value of ^0 - 4A 
is to be expected in all local "bubbles", as it 
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should be derived from the non-minimal cou- 
phng fiiR) itself: a full account of the dy- 
namics of the model, obtained by allowing the 
curvature to vary, should lead to a relaxation 
towards a constant value Rq arising from the 
modified field Eq. (2), if the composite efi'ect 
of a non-trivial curvature term and cou- 
pling f2(R) exhibits a suitable minimum, thus 
acting as an attractive potential. 

This minimum Roip) would naturally de- 
pend on the underlying density, suggesting a 
similarity with the so-called chameleon model 
(Khoury & Weltman 2004) — indeed the full 
Lagrangian /i (R) - f2{R)p may be suggestively 
compared with the usual chameleon effective 
potential Vefji4>) - y(4>) + pexp{(plM). At 
early times, local density and curvature so- 
lutions could be widely difl'erent from each 
other and from the constant curvature and spe- 
cific matter distribution assumed here: how- 
ever, extensive neighbour to neighbour inter- 
action would favour similar density profiles, 
yielding a universal minimum Rq. This would 
lead to a progressive decrease of the decelera- 
tion parameter towards q - -I, signaling the 
onset of a CC (Gong & Wang 2007). 

8. Post-inflationary reheating 

In this section one discusses the possibil- 
ity of using a non-minimal coupling to drive 
the reheating of the Universe after the in- 
flationary stage occurring in its early times 
(Bertolami et al. 2011). One does not at- 
tempt to drive the dynamics of inflation 
through the efi'ect of fiiR), but instead as- 
sumes the well-known Starobinsky inflation 
model (Starobinsky 1980) — which resorts to 
a quadratic curvature term 



m) = R- 



6M2 



(29) 



with M ~ 3 X IQ-^Mp. 

A non-minimal coupling is fundamental 
in the so-called preheating mechanism: the 
reheating of the ultracold post-inflationary 
Universe due to the explosive production of 
particles, occurring when the dynamics of a 
quantum scalar field x endowed with a vari- 
able mass term of the form m^^^ = rrp- + ^R 



experience parametric resonance (Tsujikawa et 
al. 1999). More evolved couplings also lead to 
preheating, as found in Refs. (Charters et al. 
2008; Lachapelle & Brandenberger2009). 

This hints that the non-minimally coupled 
action (1) may generaUze the preheating sce- 
nario. Given the form of the variable mass term 
Meff, one assumes a Unear coupling 



/2(^)=1+2^A . 



(30) 



Since the curvature is coupled to matter and ra- 
diation, besides the scalar field one must en- 
sure that the cosmological dynamics are driven 
by the eff'ect of the quadratic curvature term 
Eq. (29) alone — i.e. that the non-minimal cou- 
pling only intervenes during preheating, not 
before. As shown in (Bertolami et al. 2011), 
this implies a perturbative regime f2{R) ~ 1 
and that 1 < ^ < lO"*, compatible with the 
weak bound ^ ^ 10^^ obtained from con- 
siderations on solar hydrostatic equilibrium 
(Bertolami e Paramos 2008). 

Decomposing the scalar field x iiito its 
Fourier modes Xk, one finds that they foUow 
the diflferential equation 



X',' 



2k \^ I2m\^ H' 



£ / R'^ HR' 

r- -6 R' 



a'l^Xk 



Xk 

to) ± nil (de- 



using the redefinition Xk 
and the new variable 2z - M{t 
pending on the sign of f). 

After slow-roll, the Hubble parameter Hit) 
experiences an oscillatory phase, so that R ~ 
(2M^/z) cos(2z): the R" term dominates and 
one rewrites the above as a Mathieu equation, 



X'; + [Ak-2qcos{2z)]Xk = Q , 
with 

2k f /2m\^ 



^1 

z 



(32) 



(33) 



Eq. (32) is the exact same form found in usual 
preheating (Tsujikawa et al. 1999): the quan- 
tum field X experiences parametric resonance 
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as the scale factor a(t) increases, with massless 
particles produced for a couphng parameter as 
low as ^ > 3, while massive particles require 
^ > 10 — well within the range 1 < ^ < 10'* 
discussed above. Thus, one finds that, as ex- 
pected, a universal non-minimal coupling may 
successfully drive the reheating of the post- 
inflationary Universe. 

9. Conclusions and Outlook 

A non-minimal coupling between geometry 
and matter covers a broad spectrum of appli- 
cations in an elegant, natural way. Hopefully, 
this hints at its relevance in the formulation of 
a fundamental quantum theory of gravity: fu- 
ture advances on the theoretical front should 
help clarify its nature. 

The distinct expressions for the non- 
minimal coupling here considered may be re- 
garded as an approximation to a more evolved 
forms for /iC/?), each vaUd in a particular 
regime: early vs. late time, central vs. long 
range, it etc.. The latter could perhaps be writ- 
ten as a Laurent series, 

2 y . (34) 

SO that one could chart yet unprobed terms 
of the above series by assessing the dynamics 
of other phenomena and environments, where 
distinct curvatures and densities are at play. ? 
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